Analytic solution for a class of turbulence problems 
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' An exact analytical method for determining the Lagrangian velocity correlation and the diffusion 

coefficient for particles moving in a stochastic velocity field is derived. It applies to divergence- 
free 2-dimensional Gaussian stochastic fields which are stationary, homogeneous and have factorized 
Eulerian correlations. 
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Test particle motion in stochastic velocity fields is a generic problem in various topics of fluid and plasma turbu- 
lence or solid state physics jj]]. The main difficulty in determining the resulting time dependent (running) diffusion 
coefficients and mean square displacements consists in calculating the Lagrangian velocity correlation function (LVC). 
This is a very complex quantity which requires the knowledge of the statistical properties of the stochastic trajectories 
determined by the random velocity field. The vast majority of existing works employ under various guises the Corrsin 
approximation The physical parameter which characterizes such process of diffusion by continuous movements 
is the Kubo number K (defined below) which measures particle's capacity of exploring the space structure of the 
stochastic velocity field before the latter changes. In the weak turbulence case K <§; 1, particle motion is of Brownian 
type and the results are well established. In the strong turbulence case (K > 1), this structure of the velocity field has 
an important influence on the LVC and on the scaling of the diffusion coefficient in K. This influence is most effective 
in the special case of 2-dimensional divergence-free stochastic velocity fields and consists in a dynamical trapping of 
i the trajectories in the structure of the field. The existing analytical methods completely fail in describing this process 
Oh| |§ an d the studies usually rely on direct numerical simulations of particle trajectories Q, on asymptotic methods 
such as the renormalization group techniques [jD or on qualitative estimates . In recent works || a rather different 
statistical approach (the decorrelation trajectory method) was proposed for determining the LVC for given Eulerian 
correlation of the velocity field. The case of collisional particles was treated in J7|. We prove here that, in the special 
case of 2-dimensional divergence-free velocity fields, under the assumptions mentioned below, this method yields the 
exact analytical expression of the Lagrangian velocity correlation valid for arbitrary value of the Kubo number. The 
assumptions concern the statistical properties of the stochastic velocity field and are rather natural for a large class 
of physical processes. It is considered to be a stationary and homogeneous Gaussian (normal) stochastic field, either 
i-H ' static or time dependent with statistically independent space and time variations such that the Eulerian correlation 
function has a factorized structure as in Eq.(Q) below. 
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Particle motion in a 2-dimensional stochastic velocity field is described by the nonlinear Langevin equation: 

^=v(x(t),i), x(0)=0 (1) 
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where x(t) represents the trajectory in Cartesian coordinates x = (xi,X2). The stochastic velocity field v(x, t) is 
divergence-free: V • v(x, t) — and thus its two components v\ an v-i can be determined from a stochastic scalar field, 
the stream function (or potential) </>(x, t), as: 

v(x,t) = V x <Xx,i)e 2 = f A _*L)0( X( t) (2) 
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where e z is the versor along the z axis. The stochastic stream function </>(x, t) is considered to be Gaussian, station- 
ary and homogeneous. Since the velocity components are derivatives of </>(x, t), they are Gaussian, stationary and 
homogeneous as well. We assume that they have zero averages: 

(<Kx,t)> = 0, (v(x,t))=0. (3) 
The Eulerian two-point correlation function (EC) of 0(x, t) is assumed to be of the form: 
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£(x, t) = (0(0, 0) 0(x, i)) = /3 2 £(x) fc(t) 



(4) 



where j3 measures the amplitude of the stream function fluctuations and () denotes the statistical average over the 
realizations of 0(x, t). £(x) is a dimensionless function having a maximum at x = 0, where its value is £ (0) = 1, and 
which tends to zero as |x| — * oo. It actually depends on the dimensionless variable x/A, where A is the correlation 
length. h(t) is a dimensionless, decreasing function of time varying from h(0) = 1 to h(oo) = 0. It depends on 
the dimensionless ratio t/r c , where r c is the correlation time. The Kubo number is defined as the ratio of the 
average distance covered by the particles during t c to A : K = Vt c /X where V — /3/A measures the amplitude of 
the fluctuating velocity. Using the definition (0) of the velocity, the two-point Eulerian correlations of the velocity 
components and the potential-velocity correlations are obtained from E(x.,t) as: 



^ n — — ~a~2^' ^22 — —jr^2~E, E12 — a _ a _ E, (5) 
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Ei,p = —E^i = ~~q^~Ei E 2< p = — E^2 = ~q^tE 



where E v ^t) = (v i (0,0)v j (x,t)) and E^ = (0(0,0) Uj(x,t)) . 

Starting from this statistical description of the stochastic stream function, we will determine the Lagrangian velocity 
correlation (LVC), defined by: 

L ii (i) = (v i (x(0),0)u i (x(t),t)). (6) 
The mean square displacement and the running diffusion coefficient are determined by this function: 

{x 2 l {t)) = 2 [ dTL u (r) (t-r), (7) 
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Di{t) - / drLuir). (8) 
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For small Kubo numbers (quasilinear regime), the results are well established: the diffusion coefficient is 
Dql — {^ 2 /t c )K 2 . At large K the time variation of the stochastic potential is slow and the trajectories can fol- 
low approximately the contour lines of (f>(x,t). This produces a trapping effect : the trajectories are confined for long 
periods in small regions. A typical trajectory shows an alternation of large displacements and trapping events. The 
latter appear when the particles are close to the maxima or minima of the stream function and consists of trajectory 
winding on almost closed small size paths. The large displacements are produced when the trajectories are at small 
absolute values of the stream function. 

The main idea in our method is to study the Langevin equation (Q) in subensembles (S) of the realizations of the 
stochastic field which are determined by given values of the stream function and of the velocity in the starting point 
of the trajectories: 

0(O,O) = 0°, v(O,0)=v°. (9) 
The LVC for the whole set of realizations can is obtained by summing up the contributions of all subensembles fl9) : 



(t) = J J #° dv° Pi (0°) Pi (v°) Lfj (t) (10) 



where Lfj (t) is the LVC in (S) and Pi (0° ) , Pi (v° ) are the Gaussian (normal) probability densities for the initial stream 
function and respectively for the initial velocity. As shown below, there are two important advantages determined by 
this procedure: (i) the LVC can be determined from one-point subensemble averages and (ii) the invariance of the 
stream function along the trajectory can be used for obtaining the average of the Lagrangian velocity and the LVC 
in(S). 

The stream function and the velocity reduced in the subensemble (S) are still Gaussian stochastic fields but non- 
stationary and non-homogeneous and with space-time dependent average values: 



2 



**(x,t) = «x,t)) s =^^+«°^, (11) 

tf(x,t) = («,(«, t)) fl = + ^ (12) 

where () s represents the average over the realizations in the subensemble (S). These averages are determined using 
conditional probability distribution; they are equal in x = and t = to the parameters determining (S): $ s '(0, 0) = 
0°, Vf(0, 0) = v® and decay to zero as x — > oo and/or t — > oo. A relation similar to (0) can easily be deduced: 

-(£• "<£:) (13) 

which shows that the average velocity in the subensemble (S) is divergence- free: V- V s (x,i) = 0. We have thus 
identified in the zero-average stochastic velocity field a set of average velocities (labeled by <p , v ) which contain 
the statistical characteristics of the velocity field (the correlation and the constraint imposed in the problem, i.e. the 
zero divergence condition). The LVC in (S) is: 

Lg.(t) = (v i (0,0)v j (x(t),t)) s = v° (« J (x(t),*)) s (14) 

and thus the problem reduces to the determination of the average Lagrangian velocity in each subensemble. 

We consider first the static case </>( x ) ( r c — ¥ K — > oo, and the EC of the stream function E(x) is independent 
of time). The stream function is an invariant of the motion (</>(x(i)) = 4>° in each realization in (S)) and: 

(0(x(t))) s = 0° (15) 

at any time. A deterministic trajectory X(£; S) can be defined in each subensemble (S) such that the average of 
the Eulerian stream function in (S) (Eq.(|TT|)) calculated along this trajectory equals the average Lagrangian stream 
function (|l5|): 

$ s (X(t;S)) = (^(t))) s ^°. (16) 



Since the Eulerian average potential ([ll]) has the value <f>° in x = 0, the trajectory X(i; S) can be determined from 
the following Hamiltonian system of equations with <f>' s (X) as Hamiltonian function: 



dX(t;S) ( d d 



dt \dX 2 dX x 



$ S mS)) (17) 



and with the initial condition X(0; S) = 0. The trajectory in each realization in (S) can be referred to this deterministic 
(realization-independent) trajectory, x(i) = X(t; S) + <5x(f). Using the definition of the velocity (j^) and expressing the 
space derivatives as derivatives with respect to the deterministic part of the trajectory in each realization and then 
averaging in the subensemble (S), the average Lagrangian velocity in (S) is obtained as: 

(v(x(£))) s = V s (X(i;S)) (18) 



where Eq.(13) was used. Thus, the average Lagrangian velocity in the subensemble (S) is just the corresponding 
Eulerian quantity calculated along the deterministic trajectory X(t;S). Since the latter is determined by solving 
Eq. (|l7|) , where the r.h.s. is the average Lagrangian velocity, it follows that X(t; S) is precisely the average trajectory 
in (S) : 

X(t;5) = (x(t)> 5 . (19) 

Similar results are obtained in the time-dependent case 0(x, t) (finite r c and K) if the space and time dependences 
are statistically independent in the sense that the EC of 0(x, t) is given by Eq. (H) . The stream function is not a true 
invariant of the motion. However the velocity is still perpendicular to V</>(x(f), t) at any moment and only the explicit 
time-dependence contributes to the variation of <f> along the trajectory: 
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Due to the factorized EC (f|) considered here, the average Eulerian stream function and velocity (fy]), ( |l2] ) can be 
written as: 

$ 5 (x,t) = $ 5 (x)/i(i), (21) 

V s (x,i) = V s (x)/i(t) (22) 

where $ s (x) and V s (x) are the corresponding quantities in the static case. We define in (S) a deterministic trajectory 
X(t; S) as the solution of the time dependent Hamiltonian system with $ S (X, t) as Hamiltonian function. Performing 
the change of variable t — > 8(t) — J * dr h(r), the time dependent Hamiltonian system reduces to Eq.(p"7|) and thus the 
trajectory X(t; S) can be written as: 

X(t;5) = X(0(i);S) (23) 

where X((9; S) is the deterministic trajectory obtained in the static case. On the other hand, the average Lagrangian 
potential corresponding to Eq.(|2~l|) can be written as (</>(x(t), t)) s — G(t)h(t) where the factor h(t) "propagates" 
unchanged from the Eulerian average to the Lagrangian one, and G(t) is the contribution of the space dependence. 
Taking the time derivative of this equation one obtains using Eq. (|2p| ) : 



- (0(x(i),t)) s = hm . 



It follows that d (0(x(t), t)) s /dt — G(t)dh/dt and thus G(t) = <p° and the average Lagrangian stream function in (S) 
is: 

(^(t),t)) s ^4> Kt). (24) 



Using Eq.(24) and (|2l|) and the definition of X(t; S), one finds that the average Eulerian stream function calculated 



along the deterministic trajectory X(t; S) equals, as in the static case, the average Lagrangian stream function: 

4> s (X(t; S)) h(t) = (0(x(t), t)) 3 = 0° h(t) (25) 
Following the same arguments as in the static case, the average Lagrangian velocity in (S) is determined as: 

(v(x(t),t)) s -V s (X(i;5)) h(t) (26) 
and the deterministic trajectory X(£; S) is the average trajectory in (S): 

X(t;S) = (x(t)) s . (27) 

We have thus obtained the subensemble averages of the Lagrangian stream function and velocity. These averages 
of random functions of random arguments appear to be equal to the average functions evaluated at the average 
argument. We note that this surprisingly simple result which is usually wrong for stochastic functions is exact for the 
special case studied here. This property is essentially due to the very strong constraint imposed by the invariance of 
the stream function along the trajectories. 

The LVC (^) and the running diffusion coefficient (||), are determined using Eqs. ([h]) and (|l4|) as: 

Lij(t) = SijV 2 F' (K9(t)) h(t), (28) 

A 2 

D(t) = — K F (K6(t)) (29) 
where F'{t) is the derivative of the function F(t) which is defined by: 
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FIG. 1. The function F(t) denned in Eq.@ 



10 



F(t) 



dpdu u 3 exp — 



u 2 (l+p 2 ) 



X(t;u,p). 



(30) 



X(r;u,p) is the component of the average trajectory (normalized by A) along v° and the dimensionless parameters 
u= |v°| /V, p = 4>° /uj3 determine the subensemble (S). It is the solution of Eq.(17). The function in Eqs.(^8|), (|29| ) 
is defined by: 



6{t) = [ h^dh 
Jo 



(31) 



We note that Eqs.(p8|)-(pl|) represent the exact solution for the diffusion problem studied here (both for static 
and time-dependent case). Two functions are involved in the expressions for the LVC and D{t). One is the time 
dependence h(t) of the EC of the stream function. This accounts for the explicit time decorrelation and remains 
unchanged when passing from Eulerian to Lagrangian quantities. The other is the function F(t) which results from 
the space dependence of the EC of the stream function (i.e. from the Lagrangian nonlinearity) . It is obtained as an 
integral over the subensembles (S) of the average displacements along the initial velocity v°. 



The trajectories obtained from Eq.(|17]) lie on closed paths (except for <jr = which correspond to a straight 
line along v°). The size of these paths is large for small |</>°| and it decreases as \<p°\ increases. At small time r, 
X(t]u,p) = ut and F(r) = r. At large time r, the trajectory X(r;u,p) turns periodically along the corresponding 
path. The period grows with the size of the path. Thus, for a given time t, the trajectories corresponding to small 
paths (large | <jP | ) rotate many times while those along large enough paths (small \<f>°\) are still opened. Consequently, 
when calculating the integrals in Eq.(|30"|), the contribution of the small paths (large \<fi°\) is progressively eliminated 
as r increases due to incoherent mixing. As r increases, smaller and smaller intervals of (jP around <fi° = effectively 
contribute to the function F(t) which consequently decays to zero as r — > oo. Thus, the function F accounts for the 
dynamical trapping of the trajectories. The characteristic features of this self-consistent trapping process observed in 
the numerical simulations are recovered in the structure of this function. It shows that only a part of the trajectories 
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(which are not yet trapped) effectively contribute to the value of the diffusion coefficient at that moment. These are 
the particles that move on the large size contour lines of the stochastic stream function which correspond to (f> = 0. 
The trapping process is evidenced at large time and determines the decay of F as F(t) ~ t~ q . The function F(t) 
obtained for £(x) = 1/(1 + x 2 /2A 2 ) is presented in Fig.l where the two regimes are observed. The value of a for this 
case is a — 0.38. The modification of £(x) determines a variation of a around this value. 

In a static stream function (K,t c — > oo, h(t) = 1) the particles move along the "frozen" contour lines of <fi(x) and 
the process is subdiffusive. This can easily be recovered in the general result ( ]28|)- (|3l"|). The average trajectories in 
the subensembles (S) are periodic functions of time and the diffusion coefficient (|29|) is D(t) = f3 F (Ft/A) . It goes to 
zero when t — ► oo as D(t) = (3 (Vt/A)~ a and the mean square displacement is (ar(£)) ~ i 1- ". 

In the time dependent case, the time variation of the stream function determines a decorrelation effect which leads 
to a diffusive process. This is reflected in the average trajectories in the subensembles (S) (determined by Eqs.(^7|), 
(p3|)) which are not anymore periodic functions of time but all of them saturate as t — > oo (possibly after performing 
many rotations around the corresponding paths). Consequently, the decay of the function F saturates at F(K 9^) 
where 6*00 is a constant of order 1 obtained as the limit t — > 00 of Eq. (|3l|) and the asymptotic diffusion coefficient is: 

A 2 

D = —KF(K9 00 ). (32) 

Tc 

In the limit of small K, the quasilinear result is recovered from Eq.(|32]) and at large K when trapping is important 
D « (A 2 /r c )iv ^V 

In conclusion, we have presented here an exact solution for the turbulent diffusion problem for a class of velocity 
fields. We have obtained analytical expressions for the LVC and D(t) which are valid for arbitrary values of the Kubo 
number and describes the complicated process of dynamic trajectory trapping in the structure of the stochastic field. 
The basic idea of the method consists of determining the LVC by means of a set of average Lagrangian velocities 
estimated in subensembles of realizations of the stochastic field which are defined taking into account the invariants 
of the motion. It can be extended, at least as a new type of approximation, to other types of stochastic velocity fields. 
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